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Centre sector transitions in QCD-like theories with dynamical quark matter are investigated. In 
the hadronic phase, these transitions still take place in the infinite volume at zero temperature limit 
despite of the explicit breaking of the centre symmetry by the matter fields. This finding is sup- 
ported by simulations of the SU(2) Yang-Mills theory with Higgs matter The phenomenological 
impact of the centre sector transitions for dense but confined matter is explained: centre dressed 
quarks acquire Bose statistics and form a so-called Fermi-Einstein condensate. This mechanism 
is further illustrated in the Schwinger model where it solves the Silver-Blaze problem. 
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1. Introduction: 

Under normal conditions, eventually the most prominent feature of hadronic matter is colour 
confinement: quarks merely act as constituents of hadrons and qualify as auxiliary fields in the 
sense that the QCD partition function is solely given in terms of states with N-sAity zero. Suppos- 
edly, the situation changes under extreme conditions, temperature and/or density, where quarks and 
gluons are liberated. This picture is corroborated at zero density and high temperatures by means 
of lattice Monte-Carlo simulations, which offer a first principle approach with controllable error 
margins. The picture is clear-cut in the heavy quark limit: the hadronic phase is characterised by 
a Wigner-Weyl realisation of the centre symmetry while the so-called gluon plasma phase at high 
temperature relates to a spontaneously breaking of this symmetry |[l], ^]. Quite recently [^, ^, it 
was argued that this picture extends to the theory with dynamical quarks as well: although centre 
symmetry is explicitly broken by the quark matter, centre sector transitions do still take place in the 
hadronic phase (see below for further details). The role of the quarks here is to merely induce 
a bias towards the trivial centre sector. Only at high temperature, the centre symmetry also breaks 
spontaneously giving rise to the quark gluon plasma phase. 

The properties of cold, but dense QCD matter is far less understood due to the infamous sign- 
problem which prevents Monte-Carlo simulations based upon Importance Sampling. A promising 
attempt to reach beyond the scope of Importance Sampling and to simulate QCD in this region of 
the phase diagram resorts to Complex Langevin dynamics ^ though further studies are needed to 
understand their convergence properties [^,^. The intuition on dense fermionic matter arises from 
the free Fermi gas model: Antiperiodic boundary conditions in (Euclidean) time direction imply 
the formation of a Fermi sphere for finite chemical potential. The spin-statistics theorem [^, ^] 
ties particles with half-integer spin to fermion statistics. Note, however, that the theorem requires 
propagating (finite mass) particles with a positive definite norm. The statistics of confined quarks 
is therefore not necessarily restricted by the spin-statistics theorem. In fact, it was shown in [ [l0| ] 
that in the confinement regime a certain "large" change of the gluon background field can be traded 
in for a change of the quark boundary conditions [[To|]. Another well known example originates 
from QCD perturbation theory: fermionic ghost fields inherit periodic boundary conditions from 
the gluon sector and evade the spin-statistics theorem since they involve negative norm states. The 
observation that quarks do not necessarily comply with anti-periodic boundary conditions, but, at 
least for an even number of colours, might acquire Bose statistics has a tremendous phenomenolog- 
ical impact: if, by virtue of the gluonic background, quarks develop a Matsubara zero mode, they 
might underdo condensation if the chemical potential equals the quark mass gap. This phenomenon 
has been called Femii-Einstein condensation [[ll], 

Considering confined quarks as auxiliary fields, a constraint to thermodynamical observables 
emerge form the sheer fact that the sole excitations are provided by hadrons: at zero temperature, 
observables are independent of the chemical potential as long as it is smaller than the hadronic 
mass gap. A widespread problem is that approximate methods and QCD models fail to comply 
with this physical fact. Accordingly, this problem has been called the Silver Blaze problem by 
Cohen in | p2| ] reminiscent of the corresponding short story by Sir Arthur Conan Doyle, in which a 
dog doing nothing provides the essential clue to solve the crime mystery. 



2 



Dense confined matter 



Kurt Langfeld 



Figure 1: Sketch of a centre sector transition from a generic gluonic lattice configuration (left) to its centre 
transform (right). 



2. Centre sector transitions 



Let us adopt the lattice notation to illustrate the impact of centre sector transitions. A generic 
gluon configuration is given in terms of the links {x) G SU (Nc) and is illustrated in figure [l], left 
panel. A particular Polyakov line is defined by 



Pix) 



N, 



tr Yl Uo{xo,x) . 

xo=l 



(2.1) 



A centre transformed configuration is generated by multiplying the time-like links of a particular 
time slice f by a centre element Z„ = z„ 1 € SU (Nc): 

Ul^"\t,x) = Z„Uo{t,x) , Vx, Uj{'\x) = U^{x) else. (2.2) 

It is easy to see that any lattice action A of pure Yang-Mills theory, which is built upon closed 
contractible loops, is invariant under the transformation (22) while the Polyakov loop is not: 



in) 



(n) 



(x) = ZnP[U^]{x) 



In 

z„ = exp J —n 



Using the spatial Polyakov loop average p, it is straightforward to associate a so-called centre sector 
to each lattice configuration: 



C{p) =n, 



n : 



aig{p) 



2nn 



mm 



(2.3) 



where 



arg{p) = (P€]0,2k] , p = \p\ exp{i(p} . 
It is straightforward to show that the centre transformation ( |2.2| ) shifts the centre sector by n: 



C{p[U^,]) = m 



Cip 



u, 



in) 



m + n . 



It is the symmetry which is spanned by the transformation (2.2) which is spontaneously broken 
in high temperature pure Yang-Mills theory [[T], ^. The phenomenon related to this breakdown is 
deconfinement: the static quark antiquark potential is no longer hnearly rising at large distances. 
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Figure 2: Right: Sector transition probability T for the pure SU(2) YM-theory and the SU(2)-Higgs theory. 
Left: Roberge-Weisz transformation undoing a gluonic centre sector transition. 



and the gluons contribute black body radiation to the thermal energy density suggesting gluon 
Uberation. 

How does this picture change for QCD-like theories, i.e., for gauge theories with dynamical 
matter which transforms under the fundamental representation of the gauge group? Let us consider 
for the moment SIJ (Nc) Yang-Mills equipped with a quark determinant arising from integrating 



over the dynamical quark fields. The transformation ( |2.2| ) is no longer a symmetry of the action 
since the quark determinant provides a bias towards the trivial centre sector n = Nc- The question 
here is whether the lattice configurations are confined to the trivial centre sector because of this 
explicit centre symmetry breaking: 



C{p[U^])=Nc, ^ Ui;'\ n^Nc 



(YM-theory with matter). 



This question was thoroughly studied in [^] for the SU(2)-Higgs gauge theory for which the Higgs 
field plays the role of the dynamical matter. To detect whether centre sector transitions do occur, we 
divide the spatial lattice universe into two parts and calculate the spatially averaged Polyakov line 
over each part. Let us call the result /?< and py. Equally well, we can associate a centre sector to 
each half of the lattice universe by the mapping C(/7<) and C{py). We then consider the probability 
T that both centre sectors are different, i.e., C{p<) ^C{py). If centre sector transitions do occur, 
C{py) and C{p<) roughly sample all sectors ranging from \ . ..Nc leaving us with T 1 — l/Nc- 
If, on the other hand, centre sector transitions are prohibited, we would have T = 0. The so-called 
"tunnelling coefficient" T is shown in figure right panel, for the pure SU(2) Yang-Mills theory as 
function of the Wilson-j8 parameter. While T is close to 1/2 at low temperatures, it rapidly drops 
for large j8 values, which correspond to the gluon plasma phase. Also shown are the findings for 
the SU(2) Higgs theory: note that T is still close to 1/2 at small j8 indicating that centre sector 
transitions do frequently occur. It is only above a certain critical value when these transitions cease 
to exist due to a spontaneous breaking of the centre symmetry besides of its explicit breaking. Also 
note that the critical j8 value for deconfinement is smaller for the Higgs theory than for the pure 
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YM-theory. This shows that the deconfinement critical temperature is smaller for the Highs theory, 
which is expected due to the matter effects. 

Let us now include dynamical quark matter q{x) which is subjected to a centre transformed 
gluon background field ?7^"' (2.2) (see figure ^ right panel). The matter gluon interaction is as- 



sumed to be of next-to-nearest-neighbour type: q{x) U^{x)q{x + jJ.). It turns out [ 10 1 that the centre 



transformation of the gluon fields can be reversed by transforming the quarks fields via 

q^"\x) = Znq{x) forX(i>t, q^"\x) = q{x) else. (2.4) 

This almost appears to be just a change of variables and therefore an invariance of the partition 



function. Note, however, that with the transformation ( |2.4| ) the boundary conditions of the quarks 
change from antiperiodic to Z,, -periodic: 

q{xol/T,x) = (-1) Zn q{xo,x) 

Let us DetAp denote the quark determinant with antiperiodic boundary conditions and let Det(„) 
represent the determinant with Z„-periodic boundary conditions. Since the partition function of 
QCD (or the QCD-like theory) is a gluonic ensemble average over all centre sectors, the partition 
function can be written as 

I m^DetAAU^] c^-SymP^,]} = I ^f/f^' (^^^Hn)[uli'''^?j expl-^YMlt/i"^^^]} ■ 

(2.5) 

Note that for an even number of colours, there is the centre element = ~1- Hence, there is a 
particular centre sector, which gives rise to periodic boundary conditions for the quark determinant 



when the gluon fields are mapped to the trivial sector by means of the transformation ( |2.4[ ). In 



zeroth order perturbation theory, i.e., for U^"' = 1, the right hand side of ( [2.5| ) given rise to a 



centre-symmetric Fermi gas model. For Nc even, it was shown in [ ]10| ] that when the chemical 
potential /i approaches the mass gap m, the baryon density p is largely determined by the centre 
sector n = Nc/2 alone: 



dE r y , {Nc even, confinement phase). (2.6) 

m e^ ^'1 — 1 



where E can be interpreted as the one-particle energy of the (modified) Fermi gas model. Notably, 
this contribution is singular for /i — )• m, and suggests an instability due to condensation quite analo- 
gously to Bose-Einstein condensation. Since the degrees of freedom which condense are fermions 
which are exposed to a non-trivial centre background field, this has been called Fermi-Einstein 
condensation (FEC). It is important to note that FEC can only occur in the confined phase: in this 
phase, transitions from the trivial centre sector n = Nc to the sector with n = Nc/2 do occur. Under 
extreme conditions, this is no longer the case: the sector transitions cease to exist, and only the triv- 
ial centre sector background is attained. In this case, the above model coincides with the standard 
Fermi gas model displaying Fermi statistics: 

f 1 

P ^ '^^ [E-^]/T _[_ I ' (high temperature deconfinement phase). (2.7) 
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The question arises whether FEC also takes place in QCD-like theories with an odd numbers of 
colours and, most importantly, in QCD. To gain first insights, the centre symmetric quark model, 
i.e., zeroth order perturbation theory, has been generalised to the SU(3) gauge group [^. The phase 
diagram can be analytically calculated in this model. It was found that, for low temperature and 
intermediate values of the chemical potential, a phase with an excess of baryon density does occur 
if the spatial volume is small enough, i.e., if the system is under pressure [||]. 



3. Lessons from the Schwinger model 



The Schwinger model [|13|], i.e., QED in two dimensions on a space-time torus, is an ideal 
testbed for tracing out new ideas since many computations can be done analytically. Here, we 
will study the phenomenology of the centre-sector transitions in this model for a non-vanishing 
chemical potential. 

The model with massless fermions was exactly solved in Hamiltonian formalism on the line 
m [0, ^ and on in [17|, 1^]. The model on the torus has been studied in [ [l9[ ] and in par- 
ticular the temperature dependence of the chiral condensate, Wilson loop correlators and Polyakov 



line correlators have been determined pOj , pl| , p2[ ]. In turns out that chiral symmetry is anoma- 
lously broken and that only states with a vanishing net baryon number appear in the spectrum. The 
latter observation is key and implies that the baryon density ought to vanish even in the case of non- 
vanishing values of the chemical potential. Non- vanishing values of the fermion chemical potential 
have been firstly considered in [23, and it was found, indeed, that the full non-perturbative par- 
tition function is independent of the chemical potential. 



The gauge potential can be written as [ |20| ] 

5i0 



2k 

—hi+diX+do^ 



(3.1) 



where the periodic functions A and (j) integrate to zero. The constant toron fields Hq i G [0, 1[ label 
the U(l) centre sectors of the model. This can be easily seen from the fact that the shift ho -^hg + a 
transforms the Polyakov line by a U(l) centre element P{x) exp{i27ra} P{x). The partition 
function factorise into a photonic part and into the centre sector average of the fermion determinant: 



Z(r,L,M) = (27r) 



det'(-A) 
det'(-A + m2) 



dho dh\ Ae.t(i$h,^) 



(3.2) 



where my is a dynamically generated photon mass, and L is the spatial extent of the torus. Assuming 
a frozen centre sector, i.e., a fixed value for /iq, would lead to a Silver Blaze problem since the 
baryon density ps is non-vanishing in this case: 



L^cc j_ r J z z* 

'^^^ nJo ^1 eP(p-t^) +z e/5(P+A') +z* 



z = exp{-27ri/io}. 



Note also that for the particular choice /zq = 1/2, we find z = — 1 implying that the one-particle 
distribution functions are of Bose type. Thus, freezing the centre sector yields the wrong physics. 
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On the other hand, the centre sector average, i.e., the integration over the toron field /jq, yields a 
result which is independent of /i ||3|], i.e.. 



and, thus, solves the Schwinger-Blaze problem in the Schwinger model. 
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